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ABSTRACT 

■^-j- \ We have performed simultaneous X-ray and optical spectroscopic observations of the 

f^S . Low Mass X-ray Binary Cygnus X-2. We have used a new data system attached to 

the Low Resolution Imaging Spectrograph instrument on Keck II to obtain spectra 
with a mean time resolution of 72.075 milliseconds, simultaneous with pointed X-ray 
observations using the PCA onboard RXTE. In this paper, we have analysed the 
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variability in both wavebands on timescales of 16 seconds. During our observations 
Cygnus X-2 covered all three branches of the Z-curve, allowing us to study how the 
changes in X-ray spectral state affect the optical emission. As the optical flux rises 
the X-ray intensity first rises on the Horizontal Branch (0 < S z < 1) but then falls 
on the Normal Branch (1 < S z < 2) and Flaring Branch (2 < S z < 3), where S z is 
a rank number characterising the position on the Z-curve. This linear increase in the 
optical flux with S z indicates the optical flux is a good predictor of the accretion rate 
(possibly normalized by its own long-term average) inferred from the Z-state S z . We 
have used this correlation to decompose the total X-ray count-rate into three distinct 
spectral components. 
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1 INTRODUCTION 



Cygnus X-2 is one of the best studied low-mass X-ray bi- 
naries (LMXBs) due to its X-ray variability and observable 
optical counterpart. Cygnus X-2 is a close binary compris- 
ing a 1.78 Mq neutron star in a 9.84 day binary orbit with a 
0.6 Mq sub-giant companion star (Orosz & Kuulkers 1999). 
Since its discovery in 1965 (Bowyer et al. 1965) there have 
been a host of observations in the optical (eg. van Paradijs 
et al. (1990)), ultraviolet (eg. Vrtilek et al. (1990)) and X-ray 
(eg. Kuulkers, van der Klis & Vaughan (1996)) wavebands. 
The optical counterpart of Cygnus X-2, V1341 Cygni, has 
a spectral type A9III and contributes about 50% of the to- 
tal visual flux (Casares, Charles & Kuulkers 1998). The re- 
processing of X-rays in Cygnus X-2 is shown by the large 
contribution of the disc to the UV continuum and emission 



line fluxes, which arise predominantly from the X-ray heated 
accretion disc (de Jong, van Paradijs & Augusteijn 1996). 

Cygnus X-2 belongs to the sub-class of LMXBs known 
as Z-sources, whose X-ray spectral states trace a Z-shaped 
pattern in X-ray colour-colour diagrams (see Hasinger & van 
der Klis 1989) . In a similar manner to the colour-colour dia- 
grams familiar to optical astronomers, whose colour indices 
are magnitude differences between two bands, X-ray colour- 
colour diagrams employ "hardness ratios", which are ratios 
of counts in a high-energy band to counts in a low-energy 
band. Typically a hard colour (hardness ratio for two high- 
energy bands) is plotted against a soft colour (hardness ratio 
for to low-energy bands). The X-ray equivalent of a colour- 
magnitude diagram is a hardness-intensity diagram, which 
plots a hardness ratio against the X-ray intensity (count s -1 ) 
in some band. The location of the source in these diagrams 
provides information about changes in the shape of the X- 
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ray spectrum and intensity, which arise for example from 
changes in the mass accretion rate onto the neutron star. 
In the case of the Z-sources, they seem to move smoothly 
along three distinct tracks, from the Horizontal Branch, the 
upper-most stroke, through the Normal Branch, the diago- 
nal stroke, and onto the Flaring Branch, the lowest stroke, 
a sequence that until recently has been interpreted as a 
monotonic increase in the mass accretion rate from below to 
above the Eddington limit (Hasinger & van der Klis 1989). 
There is however mounting evidence that this simple rela- 
tionship is not the whole story and that a further parameter 
is needed to fully describe the X-ray behaviour. This has 
been demonstrated most noticeably with the "parallel lines" 
phenomena, where the relationship between the X-ray lumi- 
nosity and QPO frequency is valid on short timescales only, 
but changes on the duration of days to weeks, indicating a 
second parameter is needed to fully characterise the X-ray 
emission of LMXBs (see van der Klis (2001) for a discussion 
of this phenomenon). 

The position of the Z in the X-ray colour-colour dia- 
gram and the shapes of the soft and hard colour-intensity 
diagrams evolve on a timescale of weeks (Kuulkers, van der 
Klis & Vaughan 1996). During this time Cygnus X-2 varies 
between three levels; a high, medium and low level, as de- 
fined by Kuulkers, van der Klis & Vaughan (1996). The low 
intensity level is rare and only observed in the binary phase 
range 0.8-0.2, indicating it may be linked to a grazing eclipse 
of the X-ray emitting region by the mass donor star. More 
normally, the system appears to move slowly between the 
medium and high intensity levels. This can be seen in the 
RXTE ASM lightcurve, a section of which is shown in Fig- 
ure 1. The source can be seen to move slowly between the 
medium level, with a count rate of ~ 30 counts/sec and the 
high level where the count rate is ~ 50 coutns/sec. However, 
large, short timescale deviations from these levels can also 
be seen. These sudden drops in the count rate are caused 
by motion along the Z-curve and are not thought to be re- 
lated to the longer timescale changes (the so-called "secular 
variations" ) . The upper edge of this envelope represents the 
count rate at the Hard Vertex, between the Horizontal and 
Normal Branches, where the count rate is highest. 

Both the colour-colour diagrams and colour-intensity 
diagrams for the medium and high intensity levels are very 
different. In the high intensity level, the hard colour inten- 
sity diagram shows linear, almost horizontal Horizontal and 
Flaring Branches, whilst in the medium level these branches 
are more curved. The Flaring Branch is also truncated dur- 
ing the High level, which has been interpreted as obscuration 
by material in the binary, whereas in the medium level, the 
Flaring Branch intensity varies smoothly with the X-ray in- 
tensity. It is thought that the truncation is most noticeable 
during the Flaring Branch as it coincides with the highest 
mass accretion rate, which causes the inner disc to swell and 
shield our view of the X-ray emitting regions (see Kuulkers, 
van der Klis & Vaughan (1996) for a more detailed descrip- 
tion). 

Cygnus X-2 exhibits Type I X-ray bursts (eg. Kuulkers, 
van der Klis & van Paradijs (1995)), confirming the pres- 
ence of a neutron star, as seen in many other neutron star 
binaries, including the Z-source GX 17+2 (Kahn & Grind- 
lay (1984), Kuulkers et al. (1997)). These bursts typically 
have rise times of 0.1-10 seconds and are thought to be 
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Figure 1. The RXTE All Sky Monitor lightcurve for Cygnus X- 
2. The vertical lines show the beginning and end of the Keck and 
RXTE observations used in our analysis. 



due to thermonuclear burning of accreted material on the 
surface of the neutron star (Smale 1998). Unlike neutron 
stars, black holes have no stellar surface and Type I X-ray 
bursts are not observed. There is also tentative evidence 
for a third periodicity in Cygnus X-2 with a period of 78 
days, although the interpretation of this is complicated by 
the changes in X-ray flux as Cygnus X-2 moves around the 
colour-colour diagram. The long-term quasi-periodicity can 
be seen in Figure 1, where the peak-to-peak separation is 
~ 30 — 40rf. Wijnands, Kuulkers & Smale (1996) associate 
this periodicity with that of a precessing accretion disk, as 
seen in Hercules X-l (Giacconi et al. 1973). 

A previous multi-wavelength campaign on Cygnus X-2 
during June and October 1988, involving X-ray (Hasinger 
et al. 1990), UV (Vrtilek et al. 1990), optical (van Paradijs 
et al. 1990) and radio (Hjellming et al. 1990) observations, 
revealed many interesting features of the correlated vari- 
ability, which are summarized schematically in Figure 9 of 
Hasinger et al. (1990). They found that there was a direct 
correlation between the UV and optical fluxes and the po- 
sition on the Z-diagram, although the sampling of the opti- 
cal and UV data (~ 100's of seconds) only allowed general 
trends to be observed. They also found a more complicated 
behaviour between the X-ray intensity and the position on 
the Z diagram. In this series of papers we will build upon 
these observations, with higher time resolution optical data 
(~ seconds), which enables us to observe the previously ob- 
served trends in much more detail. 

In this paper, the first in the series, we shall discuss the 
observations of the X-ray and optical spectral variability of 
Cygnus X-2 in terms of broadband variability on timescales 
of 10 2-5 seconds. We will also go into some depth into the 
data reduction tasks performed, which will be abridged in 
subsequent papers. In subsequent papers, we will discuss the 
variability of narrow band features in the optical spectra 
and how this is affected by the X-ray spectral state. We 
will also discuss a simultaneous timing analysis of the data 
and the nature of the rapid (~ seconds) correlated X-ray 
and optical variability in terms of the time delayed optical 
emission relative to the driving X-ray variability. 
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2 DATA 

Simultaneous Keck and RXTE observations of Cygnus X-2 
were taken during UT 1998 July 2-6 as part of our campaign 
to search for correlated optical and X-ray variability in X-ray 
binaries. In this paper, we will show the optical data from 
the first and last nights of our campaign, where atmospheric 
conditions were most favourable. An analysis of the entire X- 
ray dataset showed that using only this subset did not reduce 
our coverage of the different spectral states of Cygnus X-2. 



2.1 Optical 

The optical data were taken using the Low Resolution Imag- 
ing Spectrograph (LRIS; Oke et al. (1995)) on the 10-m Keck 
II telescope on Mauna Kea, Hawaii. A summary of the ob- 
servations is given in Table 1. The LRIS was used with a 5.2 
arcsecond slit masked with aluminized mylar tape to form 
a square aperture. The 300/5000 grating used has a mean 
dispersion of 2.55 A/pixel in the range 3600A - 9200A. 

We used a novel data acquisition system, which reads 
out the CCD continuously, to obtain more than 220,000 
spectra of Cygnus X-2, in the form of continuous byte 
streams lasting typically 50 minutes. The mean integration 
time was measured to be 72.075ms and there is no dead-time 
between individual spectra. 

In addition to the 2048-pixel spectrum, 25-pixel under- 
scan and 75-pixel over-scan regions were used to measure the 
CCD bias level, which we subtracted from each spectrum. 
The noise for a given pixel was calculated using a readout 
noise of 6.3 e~ and a gain of 4.7 e~/ADU. Cosmic rays were 
rejected with a threshold of 10 sigma from the de-biased 
frames. A master Hatfield image was created by finding the 
median of 700 individual Hatfield spectra. This image showed 
no deviations above 0.3 % in all but 3 pixels. It was decided 
that it was therefore not necessary to Hatfield individual 
spectra. Calibration arc spectra and spectra of the sky in 
the region of the object were taken at regular intervals. 

The background spectrum, which accounts for ~ 0.9% 
of the total flux, conains two components; one from scattered 
light within the spectrograph and a second from the sky 
brightness. The scattered light level was measured from the 
counts in the pixels below 3400A, where both the GRISM 
and atmospheric transmission are negligable. This level was 
removed from all spectra and then the sky spectrum was 
subtracted using the best fit to spectra taken of a blank sky 
field. These sky spectra were taken at the beginning and 
end of the exposure, so that long timescale variations could 
be detected. The mean and variable components of these 
spectra were found by creating a lightcurve for each pixel 
and extracting the mean and gradient of this lightcurve. In 
order to supress the pixel-to-pixel variations, the mean and 
gradient were then averaged in wavelength, with a running 
median filter of width 101 pixels. The gradient in time of the 
data was found to be consistent with zero, indicating that 
the spectrum was constant for the duration of the exposure 
(~ 50 minutes). This spectrum is more than a simple sky 
spectrum and contains both the spectrum of the sky back- 
gound and an additional component thought to be due to 
scattered light within the spectrograph. 

The wavelength calibration was accomplished by fitting 
a second-order polynomial to 7 lines in a median spectrum 
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Figure 2. The average X-ray (top panel) and optical (bottom 
panel) spectra from observation 'G'. The X-ray spectrum is plot- 
ted with a logaritmic scale. 



of exposures of Hg and Ar arc lamps. Arc spectra from the 
beginning and end of the run were used to take into ac- 
count any drifts in the wavelength scale. The wavelength 
calibration was applied using the MOLLY spectral analysis 
package. 

The individual spectra of Cygnus X-2 were flux cali- 
brated using exposures taken of the standard star Feige 67 
(Oke 1990). We fitted a low-order polynomial to the median 
spectrum of the observations of this standard star. The re- 
sulting fit was used as a flux calibration and again applied 
using MOLLY. The average spectrum from observation 'C 
is shown in Figure 2. Several emission and aborption fea- 
tures can be cleary seen and will be described in detail in 
Paper II. 

2.1.1 Time calibration 

Due to the nature of the continuous stream of data, it was 
impossible to attach accurate timemarks to each individ- 
ual spectrum. Individual time marks were placed after every 
other spectrum using the computer clock. In order to estab- 
lish an absolute time reference, secondary timestamps were 
obtained on many occasions during the 5-night observing 
run. To create these timestamps, an incandescent lamp that 
illuminated the CCD was modulated in synch with WWVH 
time signals. From these timestamps we found an ephemeris 
for the observation and the mean exposure time. These were 
checked against the computer's time marks and found to be 
in excellent agreement. 

The input times were corrected to the solar system 
barycentre, using the subroutines supplied by the Interna- 
tional Astronomical Union's (IAU) working group on "Stan- 
dards Of Fundamental Astronomy" (SOFA*). Finally, these 
times were converted to seconds since JD 2450996.0, in order 
to make the numbers more manageable. 

In order to increase the signal-to-noise of the individual 

* http://www.iau-sofa.rl.ac.uk/ 
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spectra and make the data-set more manageable, the spectra 
were binned in time by a factor of 14 to give 1646 binned 
spectra with a time resolution of 1.01 seconds. 

2.2 X-ray 

The X-ray data taken with the PCA onboard the RXTE 
satellite covered 42 ksec of exposure time during UT 1998 
July 2-6. We have used a subset of this data, which amounts 
to 20.3 ksec and represents the times of simultaneous cov- 
erage during favorable atmospheric conditions. During this 
time all 5 PCUs were switched on. The data were anal- 
ysed using the ftools v5.1 software suite. The Standard- 2 
mode data were used with a time resolution of 16 seconds to 
create colour-colour and hardness-intensity diagrams, and 
to study the simultaneous variability. The lightcurves were 
extracted, using FTOOLS package SAEXTRCT, in 4 energy 
ranges, using the Standard-2 channels 1-5, 6-12, 13-21 and 
22-39. This corresponds to 1.94 - 3.72 keV, 3.72 - 6.22 keV, 
6.22 - 9.46 keV and 9.46 - 16.02 keV. The times were cor- 
rected to the solar system barycentre using the FTOOLS pack- 
age faxbary. The minimum count rate of 4000 count s _1 
implies that background subtraction was not necessary. The 
deadtime is estimated to be < 2%, which is smaller than the 
variability we are studying and therefore a deadtime correc- 
tion was not applied to the X-ray lightcurves. The average 
spectrum during observations 'G' is shown in Figure 2. 

The broadband X-ray spectral variations of Z-sources 
are most clearly shown in an X-ray colour-colour diagram. 
These are shown for Cygnus X-2 in, e.g., Wijnands & van 
der Klis (2001). The colour-colour and hardness intensity 
diagrams presented by them contain data from the 5 nights 
of observations. In this paper we present a subset of this 
data for which we have simultaneous optical observations. 
The colour-colour diagram from our RXTE observation is 
shown in Figure 3. The soft colour is defined as the intensity 
in the range 3.72 - 6.22 keV, divided by the intensity in the 
range 1.94 - 3.72 keV. The hard colour is similarly defined 
as the intensity in the range 9.46 - 16.02 keV divided by the 
intensity in the range 6.22 - 9.46 keV. 



3 THE BROADBAND X-RAY SPECTRAL 
CHANGES 

Cygnus X-2 can clearly be seen in Figure 3 to trace out 
all three branches of the Z-curve, beginning on the Hori- 
zontal Branch during the first night, then moving onto the 
lower Normal Branch and the Flaring Branch during the last 
night. If motion along the 'Z' is due to changes in the mass 
accretion rate (see Hasinger et al. (1990)), we would infer 
that M increases from night to night. 

The hardness intensity diagrams, which plot the X- 
ray colour changes as a function of intensity show that 
the hard colour varies smoothly with intensity whilst on a 
given branch, as expected from previous observations of the 
source. The transition between the Normal Branch and the 
Flaring Branch, referred to as the soft vertex, is clearly seen, 
but the Horizontal Branch to Normal Branch transition re- 
gion, referred to as the hard vertex, is not covered in any of 
our observations. 

The soft colour-intensity diagram also shows a linear 



trend on the Horizontal and Normal Branches. This is in- 
terpreted as changes in the X-ray production mechanism 
as a function of the increasing mass accretion rate. The 
soft colour-intensity diagram deviates from the expected lin- 
ear relationship (see eg. Kuulkers, van der Klis & Vaughan 
(1996)) whilst on the Flaring Branch. The colour-intensity 
diagrams are characteristic of the high intensity level colour- 
intensity diagram described by Kuulkers, van der Klis & 
Vaughan (1996), exhibiting an almost horizontal FB in the 
hard colour-intensity diagram. 



4 THE SIMULTANEOUS X-RAY AND 

OPTICAL VARIABILITY IN CYGNUS X-2 

4.1 General features of Night 1, UT 01-02/07/98 

The simultaneous X-ray and optical observations of 
Cygnus X-2 from Night 1 of our observations are shown 
in Figure 4. The position in the X-ray colour-colour dia- 
gram (the pluses in Figure 3) show that Cygnus X-2 was on 
the Horizontal Branch throughout the night. It began the 
night near the Hard Vertex, but by the end of the night had 
moved along the Horizontal Branch away from the Hard 
Vertex, which coincides with the lowest X-ray and optical 
fluxes. During these observations there is a general correla- 
tion between the X-ray and optical flux. In contrast to the 
Normal and Flaring Branch lightcurves shown in Figure 5 
and described in the next section, no rapid (10 0-2 seconds) 
correlated variability is seen on the Horizontal Branch. Even 
the X-ray burst observed around t = 93, 000 sees shows no 
optical counterpart. (Correlated, time-delayed variability on 
these timescales will be analysed in a subsequent paper in 
this series.) 



4.2 General features of Night 5, UT 05-06/07/98 

At the beginning of Night 5, Cygnus X-2 was on the Flar- 
ing Branch. During the night it moved off the Flaring and 
onto the Normal Branch, finishing approximately half-way 
between the Soft and Hard Vertices. The overall X-ray inten- 
sity is now anti-correlated with the optical flux, as shown in 
Figure 5. Although this might seem surprising, it is in fact 
the same behaviour as reported by Hasinger et al. (1990) 
between the X-ray intensity and UV flux. Their interpreta- 
tion for this was that the reprocessing material in the disk, 
which gives rise to much of the observed optical/UV emis- 
sion, intercepts X-rays from a number of different lines of 
sight, whereas our observations are limited to the line of 
sight from the neutron star to the detector. As these lines of 
sight may have quite different fluxes, averaging over many 
of them gives a better indication of the bolometric X-ray 
flux and hence M. They note that the observed X-ray flux 
is sensitive to both the changing spectrum of the source and 
any anisotropy in the X-ray emission. Perhaps the existence 
of anisotropic emission could also explain the lack of an op- 
tical counterpart to the X-ray burst, as the burst is seen by 
the observer, but not the region of the binary responsible for 
reprocessing the X-ray photons into optical/UV photons. 
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Filename Observation Night XTE window UT at start Elapsed time Number of unbinncd 

(sees) spectra 



K0010 


A 


1 


1 


10:48, 02/07/98 


1654 


23268 


K0013 


B 


1 


2 


13:15, 02/07/98 


3930 


58086 


K4005 


C 


5 


1 


09:08, 06/07/98 


1460 


20230 


K4007 


D,E 


5 


2,3 


10:05, 06/07/98 


3890 


46942 


K4009 


F 


5 


4 


11:42, 06/07/98 


3746 


51926 


K4011 


G 


5 


5 


13:17, 06/07/98 


1710 


23842 



Table 1. Summary of the Keck II optical observations. The binning factor of 14 used in the lightcurves gives an integration time of 1.01 
seconds for the resulting spectra. 
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Figure 3. Top left, X-ray colour-colour diagram for Cygnus X-2, based on the RXTE PCA data between UT 01/07/98 - 06/07/98. The 
data are taken from the Standard-2 mode, with a time resolution of 16 seconds. The soft colour is defined as the ratio of the intensities 
in the ranges 1.94 - 3.72 keV and 3.72 - 6.22 keV. The hard colour is the ratio of the intensities in the ranges 6.22 - 9.46 keV and 9.46 - 
16.02 keV (an increase in colour means the spectrum hardens). Bottom left, the hardness intensity diagram for the hard colour and the 
total intensity. Top right, the hardness intensity diagram for the soft colour and the total intensity. The Horizontal, Normal and Flaring 
Branches are indicated by pluses, dots and open circles respectively. The vertical line shows the X-ray intensity at the Normal/Flaring 
Branch vertex. 
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Figure 4. Top panels, the total X-ray intensity during the 2 RXTE visits (labeled A,B) from night 1. Bottom panels, the optical 
continuum flux in the range 5000 - 5800 A. The time is seconds elapsed since JD 2450096.0, the beginning of the night. The time 
resolution of the lightcurves is 16 seconds. The relative sizes of the left and right panels has been chosen to keep the scale the same for 
both lightcurves. The X-ray burst can clearly be seen at t = 93, 000 sees, with no optical counterpart. 



5 THE CORRELATIONS BETWEEN THE 
X-RAY AND OPTICAL DATA 

In Figure 6 we plot the X-ray intensity in the soft (top panel) 
and hard (bottom panel) bands as a function of the optical 
continuum flux. A striking feature is the significant range 
in X-ray intensity Ix that occurs at each optical flux (the 
error bars are typically smaller than the markers in the plot). 
This indicates there is not a simple one-to-one relationship 
between the X-ray and optical fluxes of Cygnus X-2. The 
range over which Ix varies, at a given optical flux (F v ), 
remains rather small on the Horizontal Branch (shown as 
pluses in Figure 6), but it then increases dramatically as 
the system progresses up the Normal and Flaring Branches 
(marked by dots and open circles respectively). It appears 
plausible from Figure 6 to define smooth upper and lower 
envelopes representing the extreme values of Ix that occur 



at each optical flux. These envelopes are indicated by the 
solid and dashed lines on Figure 6. 

On the basis of the optical and X-ray variability ex- 
hibited in Figure 6, we investigate here a decomposition of 
the X-ray intensity Ix into three spectral components, 7f, 
7b and I a, where the subscripts stand for Flaring, Baseline 
and Accretion respectively (see Section 5.1). This decompo- 
sition can be expressed mathematically as 



Ix{F v ) = a(t) I ¥ {F V ) + I B (F„) + I A (F V ), 



(1) 



where a(t) is a coefficient to describe the variability in 7f, 
with < a(t) < 1. The three components are discussed be- 
low and shown schematically in the lower panel of Figure 6. 
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Figure 5. Top panel, the total X-ray intensity during the 5 RXTE visits (labeled C-G) from night 5. Bottom panel, the optical continuum 
flux in the range 5000 - 5800 A. The time is seconds elapsed since JD 2450996.0, the beginning of the first night of the run. The time 
resolution of the lightcurves is 16-seconds. 



5.1 The three components 

The first component, If, with F for "flaring", is associated 
with the rapid variation between the bounds of the envelope. 
The variability of this component is mimicked in the optical 
lightcurve, as can be seen in the correlated flares and dips 
during Observations C and D, shown in Figure 5. In other 
words, on the Normal and Flaring Branches, the area be- 
tween the full and dashed lines in Figure 6 is filled in by 
the source performing diagonal (correlated) strokes on the 
timescale of 10 1-2 seconds. 

A second interesting feature in Figure 6 is the lower en- 
velope or "baseline" in the range of Ix when F v > 4.8 mjy. 
This baseline in X-ray intensity is visited repeatedly by the 
system during observations C and D in Figure 5, where it 
appears to represent a level below which Ix never goes. We 
notice in Figure 6 that this floor in Ix decreases as the op- 
tical flux rises. Furthermore, if we extrapolate that slope to 



lower optical flux it appears to predict the lowest Ix we 
observed on the Horizontal Branch, at S z = 0. We therefore 
draw a line across the diagram, and identify the X-ray in- 
tensity of this baseline component as Ib, B for "baseline", 
where Ib decreases linearly as the optical flux rises. 

The third component is simply the difference between 
the baseline intensity 7b and the lower envelope of Ix at 
each optical F v . This component first increases along the 
Horizontal Branch, and then decreases along the Normal 
and Flaring Branches. We refer to the X-ray intensity of 
this component that varies smoothly with the optical flux 
as I a, "A" for accretion, since we suggest that it is related 
to the accreted material (see discussion). 

Our decomposition is somewhat arbitrary. For exam- 
ple, we could combine Ib and I a into a single component 
with a more complex dependence on optical flux. However, 
as described above, we believe that the morphology of the 
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Figure 6. Top, a plot of the logarithm of the Soft X-ray in- 
tensity (1.94-6.22 keV) versus optical continuum flux. Bottom, 
the same plot for the Hard X-ray intensity (6.22 - 16.02 keV). The 
solid, dotted and dot-dashed lines represent the limits of the three 
components described in the text. The Normal, Horizontal and 
Flaring Branch points are represented by pluses, dots and open 
circles respectively. 



variability can cause spurious correlations, as is the case 
here. 

The variability outside of the envelope at low optical 
fluxes F v — 2.1 — 2.2mJy is not understood and more data 
will be needed in order to determine the cause of the non- 
correlated variability. 

5.2 The spectra of the components 

Comparing the upper and lower panels of Figure 6, show- 
ing the relationship between the X-ray and optical fluxes for 
the soft and hard X-ray bands respectively, we see how the 
3 components depend on the X-ray band. The most strik- 
ing feature of this is how the hard X-ray intensity (lower 
panel) reaches its maximum at lower F v than the soft X-ray 
intensity. 

In order to determine the spectrum of each component 
and how it varies with the optical flux, we have sub-divided 
the X-ray band into a total of four bands. These bands are 
the same as those chosen to create the hardness ratios shown 
in Figure 3, namely 1.94 - 3.72 keV, 3.72 - 6.22 keV, 6.22 - 
9.46 keV and 9.46 - 16.02 keV. The results of this decom- 
position are shown in Figure 7. Ib (Top Panel) is seen to 
have a roughly constant soft X-ray flux, while the Hard X- 
ray flux drops steadily, leading to an overall softening of the 
spectrum. This is in contrast to 7f (Bottom Panel), where 
the soft flux varies most. 7a shows the least spectral vari- 
ability, showing only a change between the transition from 
the Horizontal Branch to the Normal Branch. 



optical and X-ray variations, especially the apparent floor 
in the Ix, provides sufficient motivation to further explore 
our proposed 3-component decomposition. 

There are several deviations outside this envelope, 
which we will describe on a case-by-case basis. The most 
obvious of these deviations occurs around F v — 4.3mJy, 
where Ix is below the lower edge of the envelope. This devi- 
ation can be seen in the X-ray lightcurve in Figure 5 around 
t — 434, 000 sees in the form of a dip in the X-ray lightcurve 
that has no optical counterpart. Such "non-correlated" dips 
have been seen previously in Cygnus X-2 (Kuulkers, van der 
Klis & Vaughan 1996) and interpreted as momentary obscu- 
ration by material along the line of sight to the neutron star. 
They are not expected to have an analogue in the optical 
lightcurve as the line of sight from the X-ray emitting region 
to the reprocessing region is different from the direct line of 
sight for the observer. Similar non-correlated dips are seen 
in the Flaring Branch observations, which result in the dips 
seen around F v = 4.2mJy, where the X-ray intensity drops 
to the lower "floor" level (described in the next section) . 

Another deviation, seen at the highest optical fluxes is 
associated with the reprocessed flares seen during Observa- 
tion C in Figure 5. While we would expect such correlated 
flares to remain within the envelope, as they do appear to 
have X-ray analogues, we have not taken into account the 
time delay between the X-ray and optical emission. This 
time delay, which will be quantified in another paper in this 
series, is most noticeable when the amplitude of the vari- 
ability is largest on timescales similar to the time resolution 
of the data, as is the case on the Flaring Branch. This rapid 



6 SUMMARY AND DISCUSSION 

We have obtained a unique dataset from Keck II and RXTE 
and used it to investigate the simultaneous X-ray and optical 
variability in Cygnus X-2. We have found that on timescales, 
~ 10 2 ~ 5 seconds, the X-ray and optical variability show a 
positive-correlation on the Horizontal Branch, changing to 
an anti-correlation on the Normal and Flaring Branches. 

However, on short timescales, ~ 10°~ 2 seconds , the op- 
tical lightcurves show correlated, yet time-delayed features 
of the X-ray lightcurves. This is as expected in a linear re- 
processing model. (The short timescale correlated variability 
will be discussed in more detail in a later paper.) The simul- 
taneous X-ray and optical fluxes show a complicated corre- 
lation. Initially, whilst on the Horizontal Branch they show 
a fairly tight correlation. Then, when Cygnus X-2 is on the 
Normal Branch, the X-ray and optical flux changes become 
anti-correlated, whilst still exhibiting a strong correlation. 
Then, around the middle of the Normal Branch, the cor- 
related variability changes again: the correlation becomes 
weaker with large changes in both the X-ray and optical 
variability. 

We have used the simultaneous X-ray intensity and 
optical flux measurements to characterise the spectrum in 
terms of the three components 7b, 7a and 7p, corresponding 
to baseline, accretion and flaring components. It is clear that 
there are different phases in the behaviour of Cygnus X-2, 
represented by the different branches of the colour-colour 
diagram. This would appear to represent the complex inter- 
dependency of more than one component of the X-ray emis- 
sion. 
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Figure 7. The spectra of the three identified components in the 
X-ray spectra. The components are 7b, 7a and 7p from top to 
bottom panel respectively. The spectra were calculated when the 
optical continuum flux equaled 2.3 (pluses, solid line), 2.8 (stars, 
dashed line), 3.7 (circles, dot-dashed line), 4.05 (crosses, dotted 
line), 4.4 (squares, dot-dot-dot-dashed line), 4.75 (triangle, solid 
line) and 5.1 (diamond, dashed line). The spectra for 7a when F v 
= 2.3 & 5.1 arc set to zero, as can be seen from Figure 6. 



(i) Component 7b: 

As we have said earlier, 7b appears to describe the lowest 
possible X-ray intensity, or "baseline" as function of the rank 
number, S z - The intensity is seen to fall to this level, but 
not below it, at both extremes of the Z-curve and possibly 
when the neutron star is momentarily obscured by interven- 
ing material. Component B smoothly softens between the 
observations on the Horizontal and Flaring Branches, as a 
result of a drop in intensity, predominantly in the hardest 
energy range. 

(ii) Component 7a: 

The spectral shape of 7a remains constant while on the 
Horizontal Branch and also, with a different slope, on the 
Normal and Flaring Branches, only changing in flux while on 
either branch. The spectrum is seen to change between the 
Horizontal and Normal/Flaring Branches, when the sense of 
the correlation changes, with a harder spectrum on the Hori- 
zontal Branch. Like 7b, this component is relatively constant 
on timescales of 10 2-3 seconds. However, when the intensity 



drops into this region during non-correlated X-ray dips it 
does vary on these timescales. 

(iii) Component If: 

If was defined such that it contains most of the 10 1-2 sec- 
onds variability, as expressed by a(t) in Equation 1. The 
instantaneous X-ray intensity appears to change randomly 
within its envelope, except during X-ray bursts and non- 
correlated dips, which we therefore do not associate with 
variations in 7f. The spectrum of this component shows a 
highly variable soft intensity, so that the total spectrum soft- 
ens with increasing S z . 

If, which represents the rapidly varying material is the 
softest of the three components and could be associated with 
material in a hot inner disc. At a given optical flux and 
hence inferred mass transfer rate, the flux in this component 
varies rapidly and appears to drive much of the rapid optical 
response, as can be seen by the correlated X-ray and optical 
flares in observations C and D. Such rapid variability could 
be associated with inhomogeneities in the accretion flow. 

I a , whose flux varies considerably, but with little change 
in the spectrum may be associated with the emission from 
the surface of the neutron star. Such variability could then 
be caused by the changing obscuration of the inner regions 
of the accretion flow by optically thick material. 

However, no complete physical interpretation of the de- 
composition is possible as yet. This is particularly evident 
when considering 7b, which represents the lowest ('Base- 
line') level is the spectrum of the source when the accretion 
rate is at its lowest and highest levels. At these times there 
will be contributions from many regions within the binary, 
leading to a complicated spectrum. This component may 
include the non-thermal emission, which has been observed 
to increase with decreasing M. Observations of similar vari- 
ability of the so-called hard-tail in Cygnus X-2 (Di Salvo 
et al. 2002) and other Z-sources would appear to support the 
existence of this variable hard component (Di Salvo et al. 
2000; Frontera et al. 1998). In addition to this, it has been 
observed that the radio emission, which is linked to the mass 
outflow from the system, possibly in the form of a jet (Fender 
& Hendry 2000), is greatest at the lowest S z (Penninx et al. 
1988) and hence optical flux, as seen in our observations. 

The limited amount of data make it impossible to as- 
sess the importance of effects such as the binary phase and 
the long-term periodicity on 7a and indeed the other two 
components. However, it is intriguing that such a simple 
decomposition leads to three distinct components that are 
suggestive of the three regions thought to affect the X-ray 
emission in neutron star X-ray binaries. 
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